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ABSTRACT 

Since  the  introduction  of  parameter  plane  techniques  by  Siljak 
in  1964,  much  effort  has  been  devoted  to  the  application  of  these 
techniques  to  various  control  system  problems. 

This  paper  concerns  the  application  of  parameter  plane  techniques 
to  the  design  of  active  filter  systems. 

Constant  bandwidth  curves  on  the  parameter  plane  are  compared  with 
frequency  response  curves  to  validate  the  constant  bandwidth  curves  as 
accurate  representations  of  active  filter  performance.   This  in  turn 
implies  that  the  constant  bandwidth  curves  can  be  used  in  the  design  of 
active  filter  systems. 

Also  presented  is  an  example  of  the  use  of  constant  zeta  and  omega 
curves  on  the  parameter  plane  to  manipulate  the  pole  locations  of 
individual  circuits  to  achieve  a  desired  overall  system  configuration 
with  the  proper  frequency  scaling. 
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1.   Introduction 

The  use  of  operational  amplifiers  in  active  filters  increases  the 
flexibility  with  which  signal  filters  can  be  designed.   Primary  reasons 
for  this  increased  design  flexibility  are  the  elimination  of  the  need 
for  inductors  and  of  the  need  for  impedance  matching.   Active  filter 
circuits  can  be  used  over  a  wide  frequency  range,  oftentimes  with  the 
added  benefits  of  less  cost,  smaller  size  and  less  weight.   2   These 
features  suggest  that  active  filters  provide  a  suitable  answer  to  a  wide 
variety  of  filtering  problems,  such  as,  aerospace  medicine  data  collec- 
tion and  the  analysis  of  vibrational  fatigue  in  aircraft  structures. 

This  paper  is  a  case  study  based  upon  an  active  filter  parallel-T 
filter  designed  by  Mr.  William  Kerwin  of  NASA's  Ames  Laboratory  and  is 
primarily  concerned  with  establishing  that  parameter  plane  techniques  can 
be  effectively  applied  to  design  of  active  filter  systems.   Figure  1  is  a 
circuit  diagram  of  the  filter  and  Figure  2  illustrates  the  observed 
experimental  filter  performance  of  that  circuit.   The  experimental  per- 
formance will  serve  as  the  fundamental  reference  for  the  evaluation  of 
the  parameter  plane  design  techniques  used  in  this  paper. 

Two  general  types  of  parameter  plane  curves  will  be  considered: 
(1)  constant  bandwidth  curves  and  (2)  constant  zeta  and  omega  curves. 

Hollister  developed  equations  for  displaying  curves  of  constant 
bandwidth  on  the  parameter  plane  and  defined  a  constant  bandwidth  curve 

•w 

A  constant   bandwidth  curve   for  G(j6J^)    =  M   is    a  curve  drawn 
upon  the  parameter   plane  which   specifies    the   relation   between 
the   parameters   necessary   if   the    transfer   function  G(s),   which 
is    a   function  of   the   parameters,    is   to  have   magnitude   M  at 
the   real   frequency  CUh' 
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Nutting  developed  computer  programs  to  plot  the  constant  zeta  and 
constant  omega  curves  on  the  parameter  plane.   4    If  these  parameter 
plane  curves  are  frequency  normalized,  then  frequency  scaling  must  be 
performed  when  cascading  several  circuits  together  to  form  a  filter 
system.   A  technique  for  proper  frequency  scaling  to  allow  the  manipula- 
tion of  system  pole  locations  will  be  demonstrated.   Proper  conversion  of 
parameter  plane  root  locations  into  correct  frequency  scaled  values  for 
the  individual  circuits  permits  the  individual  circuit  transfer  functions 
to  be  multiplied  together  to  obtain  the  system  transfer  function  which 
represents  the  physical  system. 

In  establishing  parameter  plane  techniques  as  applicable  to  active 
filter  design,  the  following  approach  will  be  made. 

First,  a  mathematical  model  for  the  filter  system  will  be  developed, 
that  is,  the  transfer  function  for  the  system  will  be  derived. 

Then,  a  digital  computer  program  for  the  CDC  1604  will  be  written 
to  evaluate  the  frequency  response  of  the  mathematical  model  of  the 
system.   This  frequency  response  will  be  compared  with  the  observed 
experimental  response.   If  the  performance  of  the  mathematical  model 
conforms  with  the  actual  experimental  response  then  the  mathematical 
model  will  be  considered  a  valid  representation  of  the  physical  system. 
Furthermore,  each  of  the  individual  circuits  which  together  comprise  the 
overall  system  will  be  considered  to  be  accurately  represented  by  their 
respective  individually  derived  mathematical  models,  i.e.,  transfer 
functions. 

The  frequency  response  of  the  transfer  functions  for  the  individual 
circuits  will  then  be  considered  accurate  enough  to  serve  as  standards 
for  the  evaluation  of  the  constant  bandwidth  curves  on  the  parameter 
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plane.   Good  correlation  between  the  constant  bandwidth  curves  and  the 
corresponding  frequency  response  curves  then  implies  that  both  of  the 
curves  accurately  describe  circuit  performance  and  that  the  parameter 
plane  can,  therefore,  be  used  as  a  design  tool. 
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2.   Mathematical  Modeling  of  the  System 

Evaluation  of  the  effectiveness  of  parameter  plane  techniques 
could  be  accomplished  in  several  ways.   One  method  would  be  the  Hard- 
ware approach.   Use  parameter  plane  techniques  to  determine  the  values 
of  selected  parameters  which  will  produce  a  given  performance,  then 
actually  build  the  system  and  adjust  the  chosen  parameters  to  the  speci- 
fied values.   Observe  the  experimental  performance  and  compare  with  the 
desired  performance  to  substantiate  the  validity  of  the  design  technique. 

A  second  method  is  simulation.   In  this  case  a  mathematical  model 
of  the  system  is  developed  and  the  performance  of  the  system  is  evaluated 
on  either  a  digital  or  analog  computer.   I^ien  facilities  are  available 
for  either  analog  or  digital  simulation  these  methods  have  the  advantage 
of  rapid  evaluation  and  excellent  flexibility. 

Even  if  the  first  method  of  hardware  analysis  is  applied,  the  devel- 
opment of  the  parameter  plane  curves  requires  an  accurate  mathematical 
model  of  the  system. 

Therefore  the  first  step  in  design  using  the  parameter  plane  is  to 
develop  the  mathematical  model  of  the  system. 
2.1.   Derivation  of  Transfer  Functions  Through  Signal  Flow  Graph  Analysis 

The  mathematical  modeling  of  the  filter  system  shown  in  Figure  1  will 
be  based  upon  the  analysis  techniques  of  the  signal  flow  graph  introduced 
by  S.  J.  Mason.   5   The  system  will  be  divided  into  the  three  circuits 
indicated  in  Figure  1.   Transfer  functions  for  the  individual  circuits 
will  be  derived  using  signal  flow  graph  analysis.   The  overall  system 
transfer  function  will  be  the  product  of  the  individual  circuit  transfer 
functions.   This  overall  system  transfer  function  will  be  evaluated  by  a 
computer  program  to  be  presented  in  Section  3. 
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2.1.1.   Signal  Flow  Graph  Analysis  of  Circuit  1. 
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Figure    3.      Circuit  Diagram   for   Circuit    1 


Figure  4.      Signal   Flow  Graph   for  Circuit    1 

The    loops   of   the   signal-flow  graph,    LR  " 
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Non-touching    loops   taken  two   at    a  time,    L£ 


(2) 


L1L3  ~    11   22s 


Forward  path  transmissions,  p^ 

PV=  K  (ein,  ex,  e2,  eQ) 

Forward  path  cof actors,  /\ , 
Graph  determinant,  A 
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Transfer  function,  T    (S) ,  for  circuit  1 
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.     .       2.1.2      Signal-flow  Graph  Analysis   of  Active  Parallel-T  with   Resistive 
Feedback. 

The   circuit   diagram  shown  below  is   an   active  parallel-T  network 

and  represents  both  circuit   2   and   circuit    3. 
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Figure   5.      Circuit    Diagram,    Active  Parallel-T   Network 
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Figure   6.      Signal    Flow  Graph,   Active  Parallel-T  Network 
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The   loops   of   the   signal-flow  graph,   Lk 
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Non-touching   loops  taken  three   at   a  time,    Lx_ 
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multiply  numerator  and  denominator  by  E ,  then  factor  (p  +  1)  from 

k  +  1 


numerator  and  denominator 


T(S)    = 


:[p*  ♦  i] 


2    '    'k,1J(2.  J*_  -    K)p   -    1    *      k   *   l   R 


k 

z 

where 

p   =   RCS 

Consideration  of    load   impedance,    Z, 

let 

E   -        1 

czs 

then 


T(S)    = 


K  [p2    +    l] 


(2.2-1) 
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Substitute    into   equation    (2.2-1) 
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where 


p   =  RCS 
CZ 


C      = 


2.2.   Frequency  Scaling  the  Transfer  Function 

Whenever  several  circuits  are  cascaded  it  becomes  necessary  to 
frequency  scale  the  complete  network  to  some  reference  frequency  f 
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prior  to  evaluating  the  overall  frequency  response. 

In  the  RC  elliptical  function  filters  under  consideration  the 
"notch  frequency"  is 

f   =    l  _  RC 


n   277^    27f 


Define 


d  -  ■  —     — ■ 


then 


fo    ^o 


1   _  „/..  4   a 


Id      =  —L—   =  a(V„   = 


RC       °    R0C0 


1   --4 >p  =  -^ 


RCS     R0C0S     p     pQ 
The  normalized  transfer  function  then  becomes 


TQ(S)  = 


[(k  +  l)c'+l](^f-)2  +  [J^L(2-K)  +  (k+l)c']ij-  +  1 


2 
Multiply  numerator  and  denominator  times  a 


K  [  P§  *   a2  ] 


TQ(S)  =  

[(k+l)c'  +  l]  pg  +  a  [-^-(2  -  K)  +  (k  +  l)C']p0  ♦  a2 

Frequency   scaling  the  transfer  function: 

For  example,  consider  filter  section  2 

R   =   50  K 

C   =   638   pf  (J     =  —4; —  =    3.135   x   104 


n        R2C2 


U)Q   =   2.0  x   104 


k   =   2.0 

638   pf 

K  =   2.762  a   =  — —  =    1.5675 
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TQ2(S) 


2.762  [p2    *   2.449ll 
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2.308p2    +   1.5649(.165)p0   +   2.4491 

I-'1 '62        n2    .    9   ZiZiQi 
2.  J08        Po    +   ^•i+i+yJ- 


p2    +   0.112po    +    1.0611 


1.1967 


hi 


+   2.4491 


p2   +  0.112    p      +   1.0611 


Frequency    scaling  the  transfer  function. 
Similarly   for   filter   section    3 


R  = 

C  = 

k  = 

C  = 


(a) 


50   K 
810   pf 

2 

230 
810 

2.074 


0J    =  — - =   2.469   x   104  rad/sec 

n  R3C3 


=   0.284 


^n      2 
(-T-P-)      =    1.5220 
^o 


a   =    1.2336 


T03(S) 


2.074 


pz    +    1.5220 
*_o 


] 


[(3)  (0.2840)    +   l]p2    +   1.234  [(-(-0. 074)    +   0852o]pQ   +   1.5220 
[po   +    l-522o] 


2.074 
1.8520 


p£   +   0.492po   +   0.822 

1.1198  I  pg   +   1.522o| 
P2,   +   0.493po   +   0.822 

Frequency   scaling   the  transfer   function. 
For   filter  section   1 
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T(S)    = 


Kl 


rilT22s2    +    ril^2 


1    -    K 


1  +  _L 


T22  T11  T21 


s  +  1 


For   this    section  the   scaling   factor,    a,    is    defined: 


a?   = 


t 


0    f „  r2 


*      Til  ^2  =  -tt 


22 


^  r  s2  -  r°s2  *  Ji 

'11    '22  Q2  2 


TQ(s)    = 


2 
Po 


't2 

•  o 


1-Kl 


^22  «ll  (21 


S   +    1 


T„ 


2 
Po 


2  2 

az  a^ 

o  o 


1     •  Kl  1  +  1 


1  22  'li  T21   . 


p      +    1 


f.CS)    = 

ol 


aoKl 


>2  +  r 


'     1    -   Ki      +  1        +  1       ] 

722  Til  r  21    J 


7l   =  R„C„  =  - — trr  =  0.50  x   10 


-4 


o  o   o 


2   x   10' 


*f..    =  R1C]_    =    (50   x   103)(103  x   10-12)    =   0.50   x    10 


-4 


7^2    =  R2C2    =    (80.5   x   103)  (1.333  x   10"9)    =    1.0733  x   10~4 


T^i  =  R2ci  =  (80*5  x  lo3^10"9)  =  °-805  x  10 


-4 
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ifo  (0.5   x   10~4)(0.5   x   10~4- 


a^   =  '°        -        (0-5*10--)(0.5xl0-)    .  _1_     .   0.4658 

°  '11  ^2        (0*5  x   10~4)  (1.0733  x   lO-^)      2.1466 


aQ   =   0.6825 


(0.4658)(2.000) 
VPo>    = 


>2  +  (Q.50xlO-4)r    1-2-000        + 1 . I Is   *  0 

[l.0733xl0"4  0.50xl0-4  0.805xlO-4J 


.4658 


(0.4658) (2.000) 


pg   +    (    -0-500      +   ,    +     0-500   )        +   Q >4658 
1.0733  0.805 


(0.4658) (2.000) 


p£   +    (1    +  0.6211   -   04658)p      +  0.4658 


Frequency   normalized   transfer   function  for   section   1 

Tol(Po>    -  ~ ^ 

P£   +   1.1553p0   +  0.4658 

where 

Po   =  Rocos 


=   0.5   x   10~4S 
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3.   Evaluation  of  the  Mathematical  Model 

Does  the  mathematical  model  developed  in  Section  2  accurately 
represent  the  physical  system?  Can  those  transfer  functions  be  used 
in  developing  the  parameter  plane  curves? 

To  answer  these  questions  the  digital  computer  program,  Program 
SYS  RESP,  written  in  FORTRAN  60,  was  developed  to  evaluate  system  fre- 
quency response  when  the  system  consists  of  several  active  circuits  in 
cascade. 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


JOB0600F,  NAKAGAWA,  G.R. 
PROGRAM  SYS  RESP 

THIS  PROGRAM  MAY  BE  USED  TO  EVALUATE  THF  SYSTEM  FREQUENCY 
RESPONSE  IN  "OTH  MAGNITUDE  RATIO  AND  IN  DB  FOR  A  SYSTEM 
CONSISTING  OF  UP  TO  FIVE  ACTIVr  STAGFS  WHFN  THE  NUMERATOR 
POLYNOMIAL  FOR  THE  TRANSFER  FUNCTION  OF  EACH  STAGE  CAN  BF 
REPRESENTED  FY 

XK*(B(1)  +  B(2)*^  +  R(3)*S**2) 

AND  THF  DENOMINATOR  POLYNOMIAL  RY 

A(l)  +  A(2)*S  +  A(3)*S**2  +  A(4)*S**3 


FORMAT(6A8) 
FORMAT (6A8 ) 
FORMAT (2141 
FORMAT (F10.41 
FORMAT (F10.4) 


THE  FOLLOWING  DATA  CARDS  WILL  HE  USED 

CARD  1    FIRST  LINE  OF  GRAPH  TITLE 

CARD  2    SECOND  LINE  OF  GRAPH  TITLE 

CARD  3    NSTAGFrNUMPTS 

CARD  4    W(0)»  INITIAL  VALUE  ^r  OMEGA 

CAR  J  5    XK(1)»GAIN  COEFFI  CI  FM.T  »5TAGE  1 

CARD  6    31  1,  J)  ,  NUMERATOR  COFFF  TCI  FNTS  , STAGE  1, 

IN  ASCENDING  ORDERS  OF  S#       FORMAT ( 4F 10 . 4  ) 
CARD  7    A( 1»J) .DENOMINATOR  COFFFI CI ENTS > STAGE  1, 

IN  ASCENDING  ORDFoS  OF  S       FORMAT ( 4F 1 0 . 4 ) 
0    XK(2),  B(2»J)  AMP  A(2,J)  RESPECTIVELY 
S  SEQUENCE  TO  ACCOMODATE  THE  REQUIRCD  NUMBFR 


CARD 
REKE, 
OF  5 


8,9,1 
T  THI 
AGEr 


DIMENSION  <(5),B(5,5),A(5,5),  T TITLE (12),  W(400), 
1  BREAL( 5,400) ,BIMAO( 5,400 ) ,  3MAG(5»400),  DB(5,40O), 
2AIMAG( 5,40  0),  AMAG(5,40^),XK<^),AD0(^00),4RFAL(^»40O), 
3DBSYS(9CQ ) 

DBSYS(0)=  0.000 

DBSYS  (1)  =  0.000 

W(0)=  0.0  00 

W(l)  =  0.00  0 

LABEL  =  4H 

READ  200,  (  ITITLE(  I  ),  1  =  1,6) 

READ  200,  (ITITLE(I),  1=7,12) 

READ  30,  NSTAGE»NUMPTS 

READ  50,  W(0) 

DO  60  J  =  1,  NUMPTS 

ABS(J)  =  1.0 
6  0  CONTINUE 

DO  13  I  =  l.NSTAGE 

READ   '  0,  XK( I ) 

READ  100,  ( B(  I ,K)  ,  K  =  l  ,4) 

READ  100,  ( A( I ,K) ,  K=  1,4) 

PRINT  130, XK(  I  )  ,(B( I  ,K)  ,K=1,3) 
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25 

13 


21 


30 
50 
100 
130 
150 
200 
31C 
320 


PRINT  1 
DC  25  J 
a  (  J  )  = 
BREALt I 
BIMAG( I 
3MAG( I , 
AREAL  ( 
A  I, MAG  (  I 
AMAG  ( I 
TEMP  = 
ABS( J) 
DB( I »J) 
DBSYS( J 
COMTINU 
DO  27  J 
PRINT  3 
PRINT  3 
CONTIIMU 
CALL  DR 

1  12»ltL 
CALL  DR 

1  12,1,L 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
END 
END 


50,  (  A(  I  ,K)  ,K=1  ,4) 

=  1,NUMPTS 
W(J-l)  +  0.01 

,J  )  =  B (  I  ,  1  )  -  B (  I ,3  )   *W< J )**2 
»J)  =  B(  I  ,2  )*W( J) 

J)  =  SQRTF( BRFAL ( I ♦ J)**?  +  BIMAG ( I » J J »*2 I 
I »J)  =  A ( I  ,  1  )  -  A(  I ,3)*  W(J)**2 
»J)  =  A (1,2)*  W(J)  -  A (1,4)  *  W(J)**3 
»J)  =  SQRTF  (AREAL  (I,J)**2  +  A  I  MAG ( I , J ) **2  ) 
XK(  I  )*BMAG(  It  J ) /AMAG(  I , J ) 
=  ABS(J)*TEMP 

=  2O.*LOG10F( TEMP) 
)  =  DBSYS( J )  +  DB(  I  ,J  ) 

-  50,150,5 
10, ( DB( I ,J) ,1=1,3) 
20,  DBSYS(J) 
E 

AW(  NUMPTS,W,DB5Y5,0,0,LAR^L, ITITLF* 0,10.,°,  1,0,2,  8, 
AST) 

AW(MUMPTS,W,ABS,0,G,LA3^L, IT  I TLE, 0,0, 0,1, 0,2, 8, 
AST) 
(214) 
(Flu. 4) 
(4F10.4) 

(1X,30HTHE  XK  AND  B  CO^FFTCIFNTS  ARF  ,/,4F10.4) 
(1X,30HTHE  A  COEFFICIENTS  ARF  ,/»4F10.4) 

(6A8  ) 

(5X,18HDB  OF  THE  STAGED   ,3F18.6) 
(5X,18HDBSYS  EQUALS        ,E18.6) 


JOB  0 

X 

3  4 

.00 

.9 

1.0 

.4 

1.19 

2.44 

1.06 

1.1 

1.5 

•  8 


600  NAKAGAWA    PARALLEL-T  ACTlV^  FILTER 

=  W(J)      Y  =  DB  OUT     SK  =  2.0 
00 


316 

00 

658 

67 

91 

11 

198 

220 

22 


1.1553 

0.000 
0.112 

0.000 
0.493 


1.000 

1.000 
1.000 

1.000 

1.000 
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The  system  frequency  response  plotted  by  Program  SYS  RESP  when 
the  previously  derived,  frequency  scaled  circuit  transfer  functions  are 
inserted  as  data  is  shown  in  Figure  7.   Since  it  can  be  shown  that  3.18 
kc  corresponds  to  a  normalized  omega  equal  to  unity,  comparison  of  this 
normalized  system  response  with  the  experimental  filter  performance 
recorded  in  Figure  1  clearly  indicates  that  the  mathematical  model  is 
valid,  that  is,  the  transfer  function  is  correct.   In  addition  these 
results  also  indicate  that  the  program  will  accurately  simulate  the  per- 
formance of  the  physical  system. 

Validity  of  the  system  mathematical  model  implies: 

(1)  that  each  of  the  circuits  which  together  comprise  the  system 
evaluated  above  can  in  turn  be  represented  by  its  respective  transfer 
function. 

(2)  that  the  performance  of  these  individual  circuit  transfer 
functions  can  now  be  used  in  turn  to  evaluate  the  constant  bandwidth 
curves  on  the  parameter  plane. 
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Figure  7.      Digital   Simulation  of   System  Frequency  Response 
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4.   Active  Filter  Constant  Bandwidth  Curves  on  the  Parameter  Plane 

The  computer  program  included  in  Appendix  A  was  written  in  FORTRAN 
63  by  Jean  Bow,  Programmer,  Computer  Facility,  Naval  Postgraduate  School 
and  modified  slightly  by  the  author  to  display  constant  bandwidth  curves 
on  the  parameter  plane.   The  program  is  based  upon  the  constant  bandwidth 
equations  developed  by  Hollister  and  utilizes  the  transfer  function  for 
the  active  parallel-T  network  derived  in  Section  2.1.2. 

The  constant  bandwitdh  curves  with  amplitude  contours  in  decibels 
are  shown  on  the  parameter  plane  plots  of  Figure  8-1  through  Figure  8-18. 
As  noted  earlier  each  parameter  plane  plot  of  the  constant  bandwidth 
curves  pertains  to  a  fixed  value  of  (Oh- 

To  obtain  a  sketch  of  the  frequency  response  of  the  circuit,  select 
the  values  of  the  two  parameters  on  the  abscissa  and  ordinate.   This 
establishes  an  operating  point  on  the  parameter  plane.   Enter  each  plot 
at  the  selected  operating  point  to  determine  magnitude,  M,  at  the  fixed 
value  of  CO,    represented  by  that  plot.   The  frequency  response  curve  can 
now  be  constructed  by  plotting  M  as  a  function  of  CJ,. 

Is  this  technique  valid,  that  is,  is  it  accurate  enough  to  be  con- 
sidered as  a  possible  active  filter  design  technique? 

The  computer  program  included  in  Appendix  B  was  developed  specific- 
ally to  answer  these  questions.   Written  in  FORTRAN  60,  Program 
FREQRESP  evaluates  the  frequency  response  of  the  parallel-T  active  fil- 
ter. 

Figure  9-1  is  the  graphical  output  for  this  program  with  SK,  the 
ratio  of  the  circuit  element  values,  set  at  2.0  and  C,  the  normalized 
capacitance,  equal  to  0.20.  The  amplifier  gain  is  varied  from  1.6  to 
5.4  to  produce  the  family  of  curves  shown,  i.e.,  each  curve  represents 
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The  following  information  applies  to  Figures  8-1 
through  8-18: 

X  Variable    "Amplifier  Gain,  K" 
X  Scale        1  inch  =  0.5 

Y  Variable    "Normalized  Capacitance,  C" 

Y  Scale        1  inch  =0.1 
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Figure  8-1.   Active  Filter  Constant  Bandwidth  Curves  on  Parameter 
Plane   k  =  2.0    CU   =  0.10 
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Figure   8-2.      Active   Filter  Constant   Bandwidth  Curves  on  Parameter 
Plane     k  =  2.0         CJ  *  0.20 
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Figure   8-3.      Active  Filter  Constant  Bandwidth  Curves  on  Parameter  Plane 
k   =   2.0         CO    =  0.30 
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Figure   8-4.      Active  Filter   Constant   Bandwidth  Curves  on  Parameter  Plane 
k  =  2.0         CaJ    =  0.40 
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Figure   8-5.      Active   Filter  Constant   Bandwidth  Curves  on  Parameter   Plane 
k  =  2.0  <^  =  0.50 
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Figure  8-6.      Active   Filter  Constant   Bandwidth  Curves   on  Parameter 
k  =  2.0         CU    =  0.60 
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Figure   8-7.      Active   Filter  Constant   Bandwidth  Curves  on  Parameter  Jlane 
k  =   2.0  CO  a  0.70 
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Figure  8-8.   Active  Filter  Constant  Bandwidth  Curves  on  Parameter  Plane 
k  =  2.0    CO    =  0.80 
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Figure   8-9.      Active  Filter  Constant   Bandwidth  Curves  on   Parameter  Plane 
k   =    2.0  CO  =  0.90 
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Figure  8-10.   Active  Filter  Constant  Bandwidth  Curves  on  Parameter  Plane 
k  =  2.0  OJ    =   0.95 


Figure  8-11.      Active   Filter   Constant    Bandwidth  Curves  on  Parameter  Plane 
k  =  2.0       CO   =   1.05 
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Figure  8-12.      Active  Filter  Constant   Bandwidth  Curves   on   Parameter  Plane 
k  =   2.0        CO  =    1.10 
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Figure   8-13.      Active   Filter  Constant  Bandwidth  Curves  on  Parameter  Plane 
k   =   2.0        Cd  =   1.20 
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Figure   8-1U.      Active  Filter   Constant   Bandwidth   Curves  on   Parameter  Plane 
k  =   2.0        U)  =   1.30 
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Figure  8-15.      Active  Filter  Constant  Bandwidth  Curves  on   Parameter  Plane 
k   =   2.0        CO   =   1.40 
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Figure  8-16.   Active  Filter  Constant  Bandwidth  Curves  on  Parameter  Plane 
k  =  2.0  CO   =   1.60 
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Figure  8-17.      Active  Filter  Constant   Bandwidth  Curves  on  Parameter  Plane 
k  -   2.0        CJ  =    1.80 
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Figure  8-18.   Active  Filter  Constant  Bandwidth  Curves  on  Parameter  Plane 
k  =  2.0      C  =  0.20 
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Figure    9-1.      Active   Filter    Frequency  Response   Curves 
k    =   2.0  C   =   0.20 

X  Variable  Normalized  Radian  Frequency 

X  Scale  1"   =   0.5 

Y  Variable  Magnitude  of  Output,   M    (db) 

Y  Scale  1"  =   5.0   db 
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Figure   9-2.      Active   Filter   Frequency   Response   Curves 
k   =   2.0  C   =   0.436 


X   Variable  Normalized    Radian  Frequency 

X   Scale  1"  =  0.5 

Y  Variable  Magnitude  of   Output,   M    (db) 

Y  Scale  1"  =   5.0   db 
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Figure    9-3.      Active   Filter  Frequency  Response  Curves 
k    =   2.0  C   =   0.284 

X   Variable  Normalized   Radian   Frequency 

X   Scale  1"   =   0.5 

Y  Variable  Magnitude  of   Output,   M    (db) 

Y  Scale  1"   =   5.0  db 
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a  specific  value  of  amplifier  gain.   At  small  values  of  oj  and  again 
for  6J>.>1.0  we  observe  the  expected  result  that  output  increases  as 
gain  increases.   But  the  distinguishing  feature  of  this  family  of  curves 
is  that  near  the  "notch  frequency"  of  1.0  a  resonance  peak  is  observed  at 
some  intermediate  value  of  amplifier  gain. 

To  evaluate  the  constant  bandwidth  curves,  two  operating  points  were 
selected.   The  first  operating  point  corresponds  to  Circuit  2  in  Figure 
1.   For  this  circuit 
SK  =  2.0 
C  =  0.436 
K  =  2.7  62 
The  frequency  response  curves  for  this  circuit  are  given  in  Figure  9-2. 
A  comparison  of  constant  bandwidth  against  the  frequency  response  values 
is  provided  in  Table  I. 

The  second  operating  point  selected  corresponds  to  Circuit  3  of 
Figure  1  for  which 
SK  =  2.0 
C  =  0.284 
K  =  2.074 
The  frequency  response  curves  for  Circuit  3  are  displayed  in  Figure  9-3 
and  tabulated  in  Table  II. 

Other  operating  points  could  well  have  been  selected,  but  further 
repetition  would  do  little  more  than  reemphasize  the  conclusions  that 
are  arrived  at  by  considering  Tables  I  and  II.   For  it  is  clearly  appar- 
ent that  the  constant  bandwidth  curves  on  the  parameter  plane  are  excel- 
lent predictors  of  circuit  performance. 

It  is  therefore  valid  to  use  the  constant  bandwidth  curves  for  the 
design  of  such  active  filters. 
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EVALUATION   OF    "CONSTANT   BANDWIDTH"   CURVES 


TABLE   I   -  Circuit   2    (C   =   0.436,    K  =   2.762) 


RADIAN 
FREQUENCY 

CO 


AMPLITUDE  FROM 
CONSTANT  BANDWIDTH  CURVES 


AMPLITUDE  FROM 
FREQUENCY  RESPONSE  CURVES 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.95 

1.05 

1.10 

1.20 

1.30 

1.40 

1.60 

1.80 


9.0  db 

9.3 
10.0 
11.2 
13.5 
19.5 
18.5 

6.0 
-4.4 
-12.1 
-14.7 
-9.8 
-5.7 
-3.7 
-2.5 
-1.2 
-0.4 


9.0  db 

9.4 
10.0 
11.3 
13.6 
20.0 
19.0 

6.0 
-4.4 
-12.5 
-14.0 
-9.0 
-5.6 
-3.6 
-2.5 
-1.1 
-0.3 
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EVALUATION  OF   "CONSTANT  BANDWIDTH"  CURVES 


TABLE    II    -   Circuit    3    (C   =   0.284,    K  =   2.074) 


RADIAN  AMPLITUDE   FROM  AMPLITUDE   FROM 

FREQUENCY  CONSTANT   BANDWIDTH  CURVES  FREQUENCY  RESPONSE  CURVES 


U)  A7  A3 


0.1  6.2   db  6.2   db 

0.2  6.5  6.5 

0.3  6.8  6.8 

0.4  7.1  7.1 

0.5  7.5  7.5 

0.6  7.8  7.6 

0.7  6.0  6.0 

0.8  1.8  1.4 


ft.- 

3 

6. 

,2   db 

6. 

5 

6. 

,8 

7. 

1 

7. 

,5 

7. 

8 

6. 

,0 

1. 

,8 

-6. 

,1 

-13. 

,7 

-15. 

,7 

-10. 

,6 

-6. 

.3 

-4. 

,3 

-3. 

,0 

-l< 

,8 

+0. 

,1 

0.9  -6.1  -6.4 

0.95  -13.7  -14.0 

1.05  -15.7  -15.0 

1.11  -10.6  -10.6 

1.20  -6.3  -6.5 

1.30  -4.3  -4.3 

1.40  -3.0  -3.1 

1.60  -1.8  -1.8 

1.80  +0.1  -1.0 
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5.   Active  Filter  Constant  Zeta  and  Omega  Curves  on  the  Parameter  Plane. 

The  computer  programs  developed  by  Hutting  were  the  basis  for  the 
constant  zeta  and  omega  curves  displayed  on  the  parameter  plane  of  Fig- 
ure 10.   The  transfer  functions  used  for  this  parameter  plane  plot  are 
derived  from  the  active  filter  sections  that  have  been  designated  circuit 
#2  and  circuit  #3. 

By  careful  selection  of  the  proper  values  of  amplifier  gain,  K,  and 
the  normalized  input  capacitance,  C,  it  should  be  possible  to  adjust  the 
root  locations  to  any  desired  value  of  damping  ratio,  zeta,  and  radian 
frequency,  omega,  displayed  on  the  parameter  plane  of  Figure  10. 
5.1  Pole  Locations  and  the  Design  of  the  Low  Pass  Filter 

The  ideal  low  pass  filter  which  is  frequency  response  optimized 
has  a  gain  vs.  frequency  characteristic  shown  in  Figure  11. 

Although  such  filter  performance  is  not  physically  realizeable, 
there  are  many  filters  that  provide  acceptable  approximations  to  the 
ideal.  Two  mathematically  designed  filters,  the  Butterworth  and  the 
Chebyshev,  will  be  discussed  here. 

The  Butterworth  filter  is  said  to  be  optimally  flat  at  zero  fre- 
quency, with  amplitude  related  to  frequency  by 


H 


2    _ 


M    2n 

1  +  # 


where 

A |         is   the   amplitude 

0JC        is   the  cut-off   frequency 

n  is   the  order  of   the  filter  transfer   function 

This    relationship   is   plotted    in   Figure   12. 

The   frequency  response  of  the  Chebyshev  filter   shown   in   Figure   13 
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Amplifier  Gain,  K 
Figure  10.   Constant  Zeta  and  Omega  Curves  on  the  Parameter  Plane 
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approximates   the    ideal    frequency   response   with  a    "constant   ripple"  re- 
sponse  in  the  pass   band,    confining  the  error  to   some   specified  value 
throughout   the  pass    band.      It    should   be   noted   that   the  Chebyshev  filter 
places   equal   emphasis   on  all   frequencies    in  the  pass   band,   whereas   the 
Butterworth   filter   emphasizes   the   near  zero    frequency   behavior. 

To   achieve  the   frequency   response   shown   in   Figure    12,    the  poles   of 
the  Butterworth   filter   transfer   function  displayed  on  the   s-plane   are 
uniformly   spaced  on  a   semicircle  of    radius    (O     as    shown    in   Figure    14. 

In   a   similar   fashion   the   frequency  response   of  the   Chebyshev 
filter    is    accomplished   by  placing   the  poles   of   the    filter   transfer 
function  on   an  ellipse   on  the  s-plane.      It   has   been  shown  that   the  proper 
positioning  of  the   poles   can  be   related  to  the  poles   of   a  Butterworth 
filter  of   the   same   order   by  reducing  the  real    component   of   the   respective 
roots   by   the   factor  tanhp      6 1,   without  modifying  the   imaginary   component. 
Thus 

°~Chebyshev   =    butterworth  (tanh  ft  ) 

where 

P  =  — -     sxnh       —=- 
n  fc. 

n  is  the  order  of  the  filter 

£.  is  related  to  the  amplitude  of  the  ripple 

with       -7—  the  minimum  amplitude  in  the 

Vi  +  ez 

pass  band  as  shown  in  Figure  13. 
5.2.   Manipulation  of  Circuit  Poles  into  Chebyshev  Configuration 

The  poles  of  the  active  filter  system  shown  in  Figure  1  will  now 
be  modified  utilizing  parameter  plane  techniques  to  form  a  Chebyshev 
configuration  of  poles. 
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Figure  11.   Ideal  Low  Pass  Filter  Frequency  Response 
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Figure   12.      Butterworth  Filter  Frequency  Response 


Figure    13.      Chebyshev  Filter   Frequency  Response 
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Locus  of  Butterworth  roots       ^-~         s 


Locus  of  Chebyshev  roots 
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Figure  14.   Relative  Position  of  Poles  for  Butterworth 
and  Chebyshev  Filters  (6tn  Order) 
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The  purpose  of  this  exercise  is  two-fold: 

(1)  To  illustrate  the  methods  required  for  proper  pole  relocation 
which  must  include  an  adjustment  of  frequency  scaling  factors  introduced 
through  the  normalization  of  the  notch  frequencies  that  was  one  of  the 
steps  in  producing  the  normalized  parameter  plane  of  Figure  10. 

(2)  To  demonstrate  that  the  constant  zeta  and  omega  curves  on  the 
parameter  plane  are  an  effective  means  of  obtaining  a  desired  performance 
from  an  active  filter  system. 

For  the  purposes  of  this  discussion  the  roots  associated  with  the 
low  pass  filter  circuit  #1  will  be  assumed  to  be  fixed  and  the  positions 
of  the  remaining  poles  will  be  manipulated  to  obtain  the  Chebyshev  con- 
figuration.  Another  assumption  is  that  the  notch  or  "zero  output"  fre- 
quency is  high  so  that  the  effects  of  the  transfer  function  zero  will  be 
ignored,  that  is,  only  the  poles  will  be  considered.   The  two  assumptions 
are  made  for  the  convenience  of  clarifying  the  following  presentation. 
No  serious  limitations  upon  the  design  of  active  filter  systems  are 
introduced  by  these  assumptions. 

Figure  15  is  an  s-plane  plot  of  the  pole  configuration  of  the  active 
filter  system.   The  numbers  in  parentheses  identify  the  circuit  number 
which  is  associated  with  that  particular  pole  pair.   Thus  the  poles 
nearest  the  negative  real  axis  designated  (1)  are  associated  with  cir- 
cuit #1  and  form  the  basis  for  the  placement  of  the  remaining  two  sets 
of  po  1  es  . 

Since  the  imaginary  component  of  this  pole  pair  is  equal  to  the 
imaginary  component  of  the  corresponding  Butterworth  filter,  the  poles 
designated  (1)  establish  the  radius,  (O    ,  of  the  reference  Butterworth 
semicircle. 
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Figure  15.   Active  Filter  System  Pole  Configuration 
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CO  (1)    =    CJo   sin   165°   =  0.3638 


CJ      =         COSn      =    1.406 
°  sin   165° 


Similarly  the  real   component   of   the  poles    designated    (1)    establish 
the  magnitude   of  the  ripple   by   specifying  the  value  of   tanh 

^"Chebyshev        _  (j-  (1) 


tanh^    = 


Lutterworth  ^0  cos    165° 


tanh^    =  -0.5776 =   0.4254 

(-0.9659) (1.406) 

Now  that  CO      and  tanhyS    have  been  established,    the  remaining  pole 
positions   can  be  determined. 

CO  (3)    =    COq   sin    135°    =   0.994 

CO  (2)   =   CO     sin   105°    =    1.358 
o 

<^(1)Chebyshev  =   tanh/^(1)Butterworth   =   tanh^   cos    165° 
=   -0.578 

CT(3)_.     .       .  0.4254   CO     cos    135°   =   -0.423 

Chebyshev  o 

0"(2)        ,_      ,_  0.4254    CO     cos    105°    =   -0.155 

Chebyshev  o 

The  values   of  the   damping  ratio   for   the  respective   roots,    zeta, 
will  be   used   on  the  parameter   plane   of   Figure   10   because  of   their  qual- 
ity of  maintaining  the   same  value   completely    independent   of   any   frequency 
scaling  or    frequency   normalization  that  must   now  be   performed    in  order 
to   take   advantage   of  the    information  contained  on  the  parameter   plane 
plot   of  Figure   10   which,    it  must    be   recalled,    has   been  normalized  to    a 
"notch"  frequency  equal   to    1.0. 
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To  calculate  the  values  for  zeta  the  following  relationship  must  be 
applied: 

lo-l 


c- 


\/cr2  ♦  w2 


Thus: 


.578 


£(1)  =  —p=— 

=       Cheb        v/0.334   +   0.132 


£(3) 


_  0.423 

Cheb   "   \/0.179    I   0.988 


£(2)Cheb 


0.155 


=   0.847 


=   0.392 


=  0.113 


v5^ 


24  +    1.843 


In   summary,   the  following  values   will   produce   a  Chebyshev  config- 
uration: 

Circuit    #1 

£(1)    =   0.847 
CO  (1)    =  0.364 
(T  (1)    =   -0.578 
Circuit    #3 

£(3)   =  0. 392 
OK3)    =   0.994 
0"(3)    =   -0.423 
Circuit    #2 

£(2)    =   0.113 
GJ<i2)    =    1.358 
CT(2)    =   -0.155 
Since   the  parameter   plane  of    Figure   10    is   normalized   to   a   "notch" 
frequency   of    1.0,    the  calculated   frequencies    listed   above  must   be   divided 
by  the   transfer   function   scaling   factors,    a,    discussed   in   Section   2.2. 
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The   scale   factors   for   circuit    #3  and   circuit    #2  were   1.2  34  and   1.566. 

Therefore  the  parameter  plane   frequencies,    CO        are: 

P 

0    9941 
CO    (3)    =        *         L      =   0.8056 
P  1.234 


,      ,    x  1.358 

To  obtain  these  parameter  plane  values  of  omega  and  zeta,  enter 
Figure  10  to  determine  the  required  values  of  normalized  input  capaci- 
tance, C,  and  amplifier  gain,  K.   The  following  parameter  values  will 
result: 

Circuit  #3 

K(3)    =   1.605  a(3)    =    1.234 

G(3)    =   0.100 
Circuit    #2 

K(2)    =   2.04  a(2)    =    1.566 

C(2)    =   0.107 

Refer  to  Appendix  C  for  algebraic  verification  that  a  Chebyshev 
configuration  has  been  formed  by  the  parameter  values  selected  above. 
5.3.   Evaluation  of  Response  of  System  with  Pole  Configuration  Derived 
Through  the  Parameter  Plane. 

What  is  the  effectiveness  of  the  parameter  plane  in  realization 
of  the  "desired  performance",  which  in  this  case  was  constant  ripple  in 
the  pass  band?   Program  SYS  RESP,  presented  in  Section  3  was  used  to 
conduct  this  evaluation. 

The  overall  system  transfer  function  required  as  input  data  to 
Program  SYS  RESP  demands  proper  frequency  scaling  of  the  individual 
circuit  transfer  functions  according  to  the  formula  derived  in  Section  2 
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m   ,    ,        K(p2    +   a2) 

TQ(p)    =  ~ — ' 5" 

(3C   +   Dp2    +   a    (3.0   -    1.5K)   +    3C   p   +   a1 

Circuit    #3 


To3(p)  ■  — ~i 


1.605(p2   +   1.522) 


1.30p*   +    l.lOp   +   1.522 


Circuit    #2 


m      ,   N  2.04(p2   +   2.449) 

T      (p)    =  Y- 

oZ  1.321p      +   0.408p   +   2.449 


The   fixed   pole   associated  with   circuit    #1  retains    its   original 
circuit   transfer  function: 

T,(p)=  -9316 


ol  p2   +   1.1553p   +    .4658 

The   results   of  the   Program  SYS   RESP  evaluation  of  the   active   filter 
system  designed  through  the   application  of   constant    zeta   and   constant 
omega  curves   on  the   parameter  plane   are    shown   in   Figure   16.      These  re- 
sults   indicate  that   the  design  objective  of  obtaining   constant   ripple 
performance  has   been   attained. 
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Radian  Frequency,  60 
Figure  16.   Digital  Simulation  of  Chebyshev  System  Frequency  Response 
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6.   Conclusions 

It  has  now  been  established  that  the  constant  bandwidth  curves  on 
the  parameter  plane  and  the  constant  zeta  and  omega  curves  on  the  para- 
meter plane  provide  valid  information  which  may  be  used  in  the  design 
of  active  filter  systems. 

The  constant  bandwidth  curves  can  provide  a  very  accurate  indica- 
tion of  filter  performance  which  is  directly  related  to  parameter  values, 
In  addition  the  constant  bandwidth  curves  provide  the  filter  designer 
with  a  good  qualitative  indication  of  parameter  values  which  should  be 
avoided  because  of  extreme  sensitivity  of  filter  performance  to  varia- 
tions in  those  specific  parameter  values.   Such  regions  of  increased 
sensitivity  are  easily  recognizeable  by  the  increased  density  of  the 
contour  lines  on  the  constant  bandwidth  plots. 

When  desired  filter  performance  can  be  realized  with  a  known  pole 
configuration,  the  constant  zeta  and  omega  curves  provide  a  very  effec- 
tive means  of  achieving  that  desired  performance  through  the  proper 
specification  of  the  various  parameter  values.   These  curves  can  also 
provide  the  user  with  information  pertinent  to  the  tuning  of  the  filter 
when  performance  is  affected  by  tolerances  of  the  circuit  elements  or 
the  aging  of  components. 

It  should  be  noted,  however,  that  the  frequency  normalized  para- 
meter plane  can  yield  meaningful  information  only  with  proper  frequency 
scaling.   Considerable  care  must  be  exercised  to  correctly  frequency 
scale  the  transfer  functions  when  several  circuits  and/or  systems  are 
cascaded  to  achieve  some  desired  result. 

Subject  to  these  conditions  parameter  plane  techniques  are  useful 
tools  in  the  design  of  active  filter  systems.   The  information  contained 
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on  various  parameter  plane  presentations  should  lead  to  much  greater 
insight  into  the  problem  of  correlating  parameter  values,  pole  locations 
and  filter  performance.   Used  in  conjunction  with  the  tremendous  flexi- 
bility and  freedom  that  prevails  in  the  relocation  of  active  filter  poles, 
these  parameter  plane  techniques  should  be  substantial  assets  in  further 
development  of  filter  design  theory. 
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JOB060oF>  NAKAGA'v*  -»G.R. 
PROGRAM  FKLC;  :  oD 
THIS  PROGRAM  WAS  DEVr LOPED 
RESPONSE  OF  THE  PARALLEL- T 
FUNCTION  HAVING  A  NUMERATOR 


TO  EVALUATE 

J\-\-    FREQUENCY 

A  C  T  T  V  -  r  I L  T  r: 

"'-  WITH  A  TPANSC 

Pn[_YMO^I  AL 

OF  THE  FOR^ 

K 


(P**2  +  1) 


C 
C 
C 

f 

c 
c 

c 
c 
c 
c 
c 


AND  A  DENOMINATOR  POLYNOMIAL  OP"  THE  FORM 

(  (SK+I  )*C+1  )*P**2  +  USK+l)*'(2-K)/SK+(sk+l  )*C)*P  +  I 

WHERE 

SK  IS  THE  RATIO  OF  ELEMENT  VALUES 
C   IS  THE  NORMALIZED  CAPACJTA.NCF 

EACH  PLOT  WILL  BE  FO^  A  FIXED  VALUE  OF  SK  AND  C,  BUT 

A  FAMILY  OF  CURVES  WILL  BE  PLOTTED  TO  REPRFSENT  VARIOUS 

VALwcS  OF  AMPLIFIER  GAIN  K. 


D I MENS 

lcJREAH 
2DB( 50u 


READ  2 
READ  2 
LA(  1  ) 
LA(2) 
LA  (  3  ) 
LA  (4) 
LA(  5) 
LA  (6) 
LA(  7  ) 
LA(8) 
LA  (  9) 
LA( lu) 
LAI  li  ) 
LA (  12  ) 
LA (  13  ) 
LA( 14) 
LA ( 15) 
LA( 16) 
SK  =  6 
C  =  0. 
MAXK  = 
GAIN  = 
W  (  1  )  = 
NUMPTS 
DO  100 
GAIN  = 


ION  W 

SOU  )  » 

)  9     AM 

0  u »  ( 

Go,  ( 

=  4HK 

=  4H< 

=  4HK 

=  4HK 

=  4HK 

=  4HK 

=  4HK 

=  4HK 

=  4HK 

=  4H 

=  4H 

=  4H 

=  4H 

=  4H 

=  4H 

=  4H 

■  0 

6  0 

16 

1.4 

0.00 

=  20 

I  = 

GAIN 


»BMAGN( 500 ) »AREAL ( 5CC 

(  5UJ  ) 

TLE (  I  )  • I =1 .fa  ) 

TLE (  I  )  ,  1=7,12  ) 


(5UU),LA(40),XK(40),ITITL!r(12)»A(5)»B(c), 

BIMAG( 500 ) ,6MAGN ( 500 ) » AREAL ( 5CC ) , AIMAG ( 50C 

AGN 

ITI 

IT  I 

1.6 

1.8 

2.0 

2.2 

2.4 

2.6 

2.8 

3.0 

3.2 

K3  •■ 

K3.i 

K3. 

K4. 

K4. 

K4. 

K4. 


,4 
6 
8 
0 
2 
4 
6 


0 
0 

li  MAXK 
+  0.2 


77 


XK( I )  =  GAIN 

B(l)  =  1.0 

B(3)  =  1.0 

A(l  )  =  1.0 

A(2)  =  (2.0-XK( I ) )*{SK+1.0)/S<  +  (SK+1.0)*C 

A(3)  =  1.0  +  (SK+1.0)*C 

DO  25  J  =  2,200 

W(J)  =  W(J-l)  +  0.01 

BREAL(J)   =  3(1)  -  B(3)*W(J)**2 

BIMAG(J)  =  B(2J*W(J) 

BMAGN(J)  =  SQRTF(BREAL( J)**2  + 
A(l)  -  A(3)*W(J)**2 
A(2)*W(.J)  -  A(4}*W(J)**3 
SQRTF( AREAL( J)**2  +  AIMAG(J 


25 

33 

35 
36 

37 


AREALt  J 
A  I  MAG ( J  ) 
AMAGN( J} 
TEMP  =  Xk 


I )*BMAGN (J 
DB(J)  =  8.685  8  9* 
CONTINUE 

IF  (I-  1) 


AL(J)**2  + 
/AMAGN( J ) 


BIMAG( J ) **2 ) 


2) 


MODE  =  1 
GO  TO  55 
IF  (MAXK 
MODE  =  2 
GO  TO  55 
MODE  =  3 
GO  TO  55 


LOGF( TEMP ) 
33,33,35 


-  I)  37,37,36 


GO  TO  55 
5  5  CALL  DRAW(NUMPTS»W»D3»MODE»0fLA( I) iITIT 

10,2  ,0,9,15,  l.LAST) 
100  CONTINUE 
20  0  FORMAT  (6A8) 

END 

END 


LF  »  •  2  »  :>  •  »  8  » 


JOB  0600  NAKAGAWA 
X=W(J)      Y=Dd  OUT 


PARALLEL-T  ACTIVE  FILTER 
SK  =  6.0      C  =  0 • 6  0 
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APPENDIX  C 
Algebraic  Verification  of  Chebyshev  Pole  Configuration 

The  parameter  plane  techniques  of  Section  5  should  yield  parameter 
values  which  will  result  in  the  desired  Chebyshev  pole  configuration. 
An  algebraic  evaluation  of  the  characteristic  equations  for  circuits  1, 
2,  and  3  will  now  be  conducted  to  confirm  that  such  a  Chebyshev  pole 
configuration  has,  in  fact,  been  attained. 

The  parameter  plane  produced  the  following  frequency  scaled  transfer 
functions  for  circuits  1,  2  and  3  respectively: 

m      ,   n  0.9316 

T   ,(P)    = 


ol  p2   +    1.1553p   +   0.4658 


2.04(p2  +  2.449) 

T_2(P)  =  5 

oz  1.321p2  +  0.408  +  2.449 


,  N       1.605(p2  +  1.522) 

T  _(p)  = 

°3       1.30p^  +  l.lOOp  +  1.522 


The  characteristic  equations  (denominator  polynomials)  for  these 
transfer  functions  are  of  the  form 


Ap2  +  Bp  +  C  =  0 


p2  ♦  £  p  +  £  .  0 
A      A 

Thus   the  natural   frequency,    6J   ,    for  the  equation   is 
n      V    A 


since 


't<»n  = 


_B_ 

A 
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c- 


»K 


2A  0n  2C 


and   the   imaginary  component  of  the   pole  becomes 


Circuit   1 


0)   (1)    =      0.4658      =      0.6825 
n 


£<». 


(1.155)(0.6825) 
2(0.4658) 


=   0.847 


CO  (1)    =    (0.6825)  (0.5331)      =      0.364 
Circuit    2 


?  449 

0)   (2)    =        — =      1.362 

n  1.321 


>         =      (0.408X1.362)        _ 

^  2(2.449) 


00(2)    =    M-.  362)  (0.9936)      =      1.354 


Circuit    3 


1   522 
OJ    (3)    =  *  =    1.082 

n  1.30 


>,   .  (1.100)(1.082) 

£(3) 2TT3213 =0-391 


&J(3)    =    (1.082X0.9204)    =   0.996 

Comparison  of  these   calculated   values   of   zeta   and  omega   with  the 
desired  Chebyshev  values   of    Section   5.2    substantiates   the   effectiveness 
of   the  constant    zeta   and  omega  curves   on  the   parameter   plane    for  the  de- 
sign task  of  obtaining  a  desired   active   filter  system  pole   configuration. 
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